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We investigate the astrophysical and cosmological implications of the recently proposed idea of a running 
gravitational coupling on macroscopic scales. We find that when applied to the rotation curves of galaxies, their 
flatness requires still the presence of dark matter. Bounds on the variation of the gravitational coupling from 
primordial nucleosynthesis, change of the period of binary pulsars, gravitational lensing and the cosmic virial 
theorem are analysed. 
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1. INTRODUCTION 

The flatness of the rotation curves of galaxies 
and the large structure of the Universe indicate 
that either the Universe is predominantly made 
up of dark matter of exotic nature, i.e. non- 
baryonic, and/or that on large scales gravity is 
distinctively different from that on solar system 
scales, where Newtonian and post-Newtonian ap- 
proximations are valid. The former possibility 
has been thoroughly investigated (see Ref. 
for a review) and is an active subject of research 
in astroparticle physics. The second possibility, 
however relevant, has drawn less attention. This 
is essentially due to the fact that until recently 
no consistent and appealing modification of New- 
tonian and post-Newtonian dynamics has been 
put foward. Many of these attempts 0, although 
consistent with observations, were most often un- 
satisfactory from the theoretical point of view. 
Actually, it has been recently shown that under 
certain fairly general conditions it is unlikely that 
relativistic gravity theories can explain the flat- 
ness of the rotation curves of galaxies ||. These 
conditions however do not exclude the class of 
generalizations of General Relativity that involve 
higher-derivatives. Quantum versions of these 
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theories were shown to exhibit asymptotic free- 
dom in the gravitational coupling [Q and one 
would expect this property to manifest itself on 
large scales. This possibility would surprisingly 
imply that quantum effects could mimic the pres- 
ence of dark matter J5j, as well as induce other 
cosmological phenomena |(|f?]]. One striking im- 
plication of these ideas is the prediction that 
the power spectrum on large scales would have 
more power than the one predicted by the Q = 1 
Cold Dark Matter (CDM) Model, in agreement 
with what is observed by IRAS ||. Furthermore, 
due to the increase in the gravitational constant 
on large scales one finds that the energy density 
fluctuations grow quicker than in usual matter 
dominated Fricdmann-Robcrtson- Walker models 
Moreover, one can explain with a scale- 
dependent G the discrepancy between determina- 
tions of the Hubble's parameter made at different 
scales, as suggested in ||, and studied in |J. 

Nevertheless, independently of the possible 
running of the gravitational constant in a higher 
derivative theory of gravity, it is worthwhile 
analysing the constraints on the scale-dependence 
of G from astrophysical and cosmological phe- 
nomena, where such an effect would be dominant. 
On the other hand, in the last few years there has 
been a revival of Brans-Dicke like theories, with 
variable gravitational coupling, that has led to 
a number of constraints on possible time varia- 
tions of G. Of course, some of the constraints 



2 



on G can be written as constraints on AG over 
scales in which a graviton took a time At to prop- 
agate. For instance, during nucleosynthesis the 
largest distance that a graviton could have tra- 
versed is the horizon distance at that time, i.e. a 
few ligh-seconds to a few light-minutes, approxi- 
mately the Earth-Moon distance. Such a distance 
is too small for quantum effects to become ap- 
preciable, as we discuss below. However, those 
effects become important at kiloparsec (kpc) dis- 
tances and therefore could be relevant for dis- 
cussing the rotation curves of galaxies. We shall 
actually show, for a particular theory (5|(|, that 
the rotation curves of spiral galaxies cannot be 
entirely explained by the running of G, so some 
amount of dark matter is required, which is still 
consistent with the upper bound on baryonic mat- 
ter coming from primordial nucleosynthesis. On 
the other hand, we could impose bounds on a 
possible variation of G from a plethora of cos- 
mological and astrophysical phenomena at large 
scales, although the lack of precise observations 
at those scales make the bounds rather weak JTo) . 
Of course, a difficulty in examining constraints 
on the variation of G is that in all gravitational 
phenomena the gravitational coupling appears in 
the factor GM, and hence we cannot distinguish 
a variation in G from the existence of some type 
of dark matter. 

2. ASYMPTOTIC FREEDOM OF THE 
GRAVITATIONAL COUPLING 



The main idea behind the results of Refs. 
is the scale depedence of the gravitational cou- 
pling. The inspiration for this comes from the 
property of asymptotic freedom exhibited by 1- 
loop higher-derivative quantum gravity models 
Since there exists no screening mechanism 
for gravity, asymptotic freedom may imply that 
quantum gravitational effects act on macroscopic 
and even on cosmological scales, a fact which has 
of course some bearing on the dark matter prob- 
lem H and on the large scale structure of the Uni- 
verse [^J . ft is in this framework that a power 
spectrum which is consistent with the observa- 
tions of IRAS [S and COBE can be obtained 



We briefly outline this proposal. Removing the 
infinities generated by quantum fluctuations and 
ensuring renormalizability of a quantum field the- 
ory requires a scale-dependent redefinition of the 
physical parameters. Furthermore, the removal 
of those infinities still leave the physical param- 
eters with some dependence on finite quantities 
whose particular values are arbitrary. These can 
be assigned by specifying the value of the phys- 
ical parameters at some momentum or length 
scale; once this is performed, variations on scale 
are accounted for by appropriate changes in the 
values of the physical parameters via the renor- 
malization group equations (RGEs). Thus, the 
equations of motion in the quantum field theory 
of gravity should be similar to the ones of the 
classical theory, but with their parameters re- 
placed by the corresponding 'improved' values, 
that are solutions of the corresponding RGEs. 
However, since gravity couples coherently to mat- 
ter and exhibits no screening mechanism, quan- 
tum fluctuations of the gravitational degrees of 
freedom contribute on all scales. One must there- 
fore include the effect of these quantum correc- 
tions into the gravitational coupling, G, promot- 
ing it into a scale-dependent quantity. One-loop 
quantum gravity models indicate that the cou- 
pling G(^ 2 /^ 2 ~ r 2 /r 2 ) is asymptotically free, 
where /i* is a reference momentum, meaning that 
G grows with scale Q. A typical solution for 
G{r 2 /r 2 ) was obtained in Ref. (5| setting the (3- 
functions of matter to vanish and integrating the 
remaining RGEs: 



G(r 2 /V 2 ) = Gi ab 5(r,r lab ) , 



(1) 



where Gi a b is the value of G measured in the lab- 
oratory at a length scale ri a b, and S(r, ri a b) is a 
growing function of r. In order for the asymp- 
totic freedom of G(/i 2 //x 2 ) to have an effect on 
for instance the dynamics of galaxies and their 
rotation curves, the function S(r, ri a b) should be 
close to one for r < 1 kpc, growing significantly 
only for r > 1 kpc. A convenient parametrization 
for S(r, riab) from the fit of Ref. |J in the kpc 
range is the following: 



6(r,r lab ) = 1.485 



1+/3 
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where (3 ~ 1/30, 7 ~ 1/10 and r = 10 kpc. 

We mention that a scale dependence of the 
gravitational constant also arises from completely 
different reasons in the so-called stochastic infla- 
tion formalism jlj] and that the scaling behaviour 
and screening of the cosmological constant was 
also discussed in the context of the quantum the- 
ory of the conformal factor in four dimensions 

Hi- 

In what follows we shall use the fit (1),(2) in 
our analysis of the rotation curves of galaxies, and 
extract a prediction for the distribution of dark 
matter. However, before we pursue this discus- 
sion let us present some of the ideas developed 
in Refs. [5-7]. As discussed above, the classi- 
cal equations have to be 'improved' by introduc- 
ing the scale dependence of the gravitational cou- 
pling. This method suggests that the presence of 
cosmological dark matter could be replaced by an 
asymptotically free gravitational coupling. As- 
suming that the Friedmann equation describing 
the evolution of a flat Universe is the improved 
one, then: 

H 2 d) = ^G(4§W , 0) 
3 a z l z 

where a = a(t) is the scale factor, H — a/a is the 
Hubble parameter, p m is the density of matter, 
£ is the comoving distance and I* is some conve- 
nient length scale. 

From Eq. (||) one sees that the present phys- 
ical density parameter, ys , is by construction 
equal to one. However, the quantity which is usu- 
ally referred to as density parameter is actually: 

87T Gp mo 
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(4) 



where H 0r is the present Hubble parameter for 
a given large scale distance, r = r*. This leads 
to £! as a growing function of scale, which is in 
agreement with observations for a constant p mo . 

Furthermore, from Eq. (3) one can clearly see 
the scale dependence of the Hubble parameter 
Moreover, as shown in Refs. (6) and 0, 
the power spectrum resulting from these consid- 
erations is similar to that of a low density Cold 
Dark Matter model with a large cosmological con- 
stant [B3. 



3. ROTATION CURVES OF GALAXIES 

Let us now turn to the discussion of the im- 
plications of the fit (^) for the rotation curves 
of galaxies. It is a quite well established observa- 
tional fact that the rotation curves of spiral galax- 
ies flatten after about 10 to 20 kpc from their 
centre, which of course is a strong dynamical ev- 
idence for the presence of dark matter and/or of 
non-Newtonian physics. The rotation velocity of 
the galaxy is given by the non-relativistic relation, 



G(r)M(r) 



(5) 



which approaches a constant value some distance 



from the centre, 



e.g. vfi 



220 km/s for the Milky 



Way. Assuming that the gravitational coupling 
is precisely Newton's constant Gn and imposing 
that the rotation velocity is constant, using the 
Virial Theorem at r = R = 500 kpc, one finds the 
standard expression for the mass distribution of 
dark matter: 



M N (r) = M N (R) 



R 



(0) 



Assuming instead a running gravitational cou- 
pling (||) , the condition that the rotation velocity 
is constant yields: 

M(r) = -. — T M N (r) . (7) 
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Equation (0) reveals after simple computation 
that the running of the gravitational coupling re- 
duces the amount of dark matter required to ex- 
plain the flatness of the rotation curves of galaxies 
by about 44%, assuming that galaxies stretch up 
to a distance of about 500 kpc. This result |l0| 
(see also Ref. [|l5| ) is a clear prediction of the de- 
pendence of the gravitational coupling with scale 
and, in particular, of the fit (^). Furthermore, 
since the possibility that the Galactic halo is en- 
tirely made up of baryonic dark matter is barely 
consistent with the nucleosynthesis bounds on the 
amount of baryons |l6| ] , the running of G is quite 
welcome since it reduces the required amount of 
baryonic dark matter in the halo. An entertain- 
ing hypothesis could be that precisely this effect 
is responsible for the reduction in the microlens- 
ing event rates across the halo in the direction of 
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the Large Magellanic Cloud |l7 18 with respect 
to those along the bulge of our galaxy, as reported 
by @. 

4. BOUNDS ON THE VARIATION OF G 
WITH SCALE 

In this section we constrain the variation of the 
gravitational coupling given by the fit (|2|) with 
bounds from primordial nucleosynthesis, binary 
pulsars and gravitational lensing and also discuss 
the effect that a scale-dependent G has on the 
peculiar velocity field (!(]] . 

4.1. Primordial nucleosynthesis 

As mentioned in the introduction, one could 
obtain bounds on the variation of the gravita- 
tional coupling from observations of the light el- 
ements' abundances in the Universe. Such obser- 
vations are in agreement with the standard pri- 
mordial nucleosynthesis scenario, but there is still 
some room for variations in the effective number 
of neutrinos, the baryon fraction of the universe 
and also in the value of the gravitational constant. 
For instance, the predicted mass fraction of pri- 
mordial i He can be parametrised, in theories with 
a variable gravitational coupling, in the following 
way m, 



Y p = 0.228 + 0.010 In mo + 0.327 log£ 



(8) 



where 7710 is the baryon to photon ratio in units 
of 10 -10 and £ is the ratio of the Hubble param- 
eter at nucleosynthesis and its present value, it- 
self proportional to the square root of the corre- 
sponding gravitational constant. In the fit (|J) it 
is assumed that the effective number of light neu- 
trinos is N u = 3 and that the neutron lifetime is 
t„ = 887 seconds. 

By running the nucleosynthesis codes for dif- 
ferent values of G, it was shown in Ref. [glj that 
a variation of AG/G = 0.2 on the values of the 
gravitational coupling was compatible with the 
observations of the primordial D, 3 He, He and 
7 Li abundances at ler level. 

This result will now be used to constrain the 
running of G in an asymptotically free theory of 
gravity. In a theory with a scale-dependent grav- 
itational constant, the maximum value of G at 



a given time is the one that corresponds to the 
physical horizon distance at that time. During 
primordial nucleosynthesis, the horizon distance 
grows from a few light-seconds to a few light- 
minutes, i.e. less than a few milliparsecs. At 
that scale we find AG/G = 0.07, see Eq. (§), 
which is much less than the allowed variation of G 
given in pl| . Therefore, primordial nucleosynthe- 
sis does not rule out the possibility of an asymp- 
totically free gravitational coupling. Of course, a 
light-second is about the distance to the Moon, 
and there are similar constraints on a variation of 
G at this scale coming from lunar laser ranging, 
AG/G < 0.6 ||. 

4.2. Binary pulsars 

The precise timing of the orbital period of bi- 
nary pulsars and, in particular, of the pulsar PSR 
1913+16, provides another way of obtaining a 
model-independent bound on the variation of the 
gravitational coupling |^| . Since the semima- 
jor axis of that system is just about a few light- 
seconds, the resulting limits on the variation of 
G can be readily compared with the ones arising 
from nucleosynthesis. The observational limits on 
the rate of change of the orbital period, mainly 
due to gravitational radiation damping, together 
with the knowledge of the relevant Keplerian and 
post-Keplerian orbiting parameters, allows one to 
obtain the following limit p3[|: 



AG 
~~G~ 



< 0.08 h' 



(9) 



where h is the value of the Hubble parameter in 
units of 100 km/s/Mpc. For h = 0.8, || one ob- 
tains cr = 0.1 which is more stringent than the nu- 
cleosynthesis bound, but is still compatible with 
the fit (§). 

4.3. Gravitational lensing 

Gravitational lensing of distant quasars by in- 
tervening galaxies may provide, under certain as- 
sumptions, yet another method of constraining, 
on large scales, the variability of the gravitational 
coupling. The four observable parameters associ- 
ated with lensing, namely, image splittings, time 
delays, relative amplification and optical depth do 
depend on G, more precisely on the product GM, 
where M is the mass of the lensing object. This 
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dependence might suggest that limits on the vari- 
ability of G could not be obtained before an in- 
dependent determination of the mass of the lens- 
ing object. However, as the actual bending angle 
is not observed directly, the relevant quantities 
are the distance of the lensing galaxy and of the 
quasar. Since these quantities are inferred from 
the redshift of those objects, they depend on their 
hand on G, on the Hubble constant, Hq, and on 
the density parameter, However, as we have 
previuosly seen, a scale-dependent gravitational 
coupling implies also a dependence on scale of H a 
and Qq, see Eqs. (3) and (4). This involved de- 
pendence on scale makes it difficult to proceed as 
in Ref. pq |, where gravitational lensing in a flat, 
homogeneous and isotropic cosmological model, 
in the context of a Brans-Dicke theory of gravity, 
was used to provide a limit on the variation of G: 



AG 
~G 



= 0.2 



(10) 



Since for this limit f2g = 1 was assumed, while 
in a scale-dependent model it is achieved via the 
running of the gravitational coupling, the bound 
( fio| ) contrains only residual variations of G that 
have not been already taken into account when 
considering the dependence on scale of Hq and 
f2o Of course, for models where the cosmological 
parameters are independent of scale, the bound 
( fl0| ) can be readily used to constrain the variabil- 
ity of G on intermediate cosmological scales. It 
is worth stressing that this method, besides be- 
ing one of the few available where this variabil- 
ity is directly constrained at intermediate cosmo- 
logical times between the present epoch and the 
nucleosynthesis era, it is probably the only one 
which can realistically provide in the near future 
even more stringent bounds on even larger scales 
by observing the lensing of light from far away 
quasars caused by objects at redshifts of order 
z > 1. 

4.4. Peculiar velocity field 

Since we expect the effects of a running G to be- 
come important at very large scales, one could try 
to explore distances of hundreds of Mpc, where 
the gravitational coupling is significantly different 
from that of our local scales. That is the realm 



of physical cosmology where of particular impor- 
tance is the study of the peculiar velocity field. A 
possible signature of the running of G would be a 
mismatch between the velocity fields and the ac- 
tual mass distribution, such that at large scales 
the same mass would pull more strongly. To be 
more specific, in an expanding universe there is a 
relation between the kinetic and gravitational po- 
tential energy of density perturbations known as 
the Layzer-Irvine equation (see eg. Ref. pq] ) that 
can be written as a relation between the mass- 
weighted mean square velocity v 2 and the mass 
autocorrelation function £(r), 



v 2 {r) = 2nG p b J 2 (r) , 



(11) 



where p\> is the mean local mass density and 
= Jo r dr £(r). The galaxy-galaxy corre- 
lation function can be parametrized by £(r) ~ 
(r/ro) -1 ' 8 with ro = bh^ 1 Mpc, while the cluster- 
cluster correlation function has the same expres- 
sion with ro = 20h~ 1 Mpc. This means that 
the velocity field (|ll|) should be proportional to 
(r/ro) ' 2 , unless the gravitational constant has 
some scale dependence. So far the relation seems 
to be satisfied, under rather large observational 
errors (see Ref. [^7j for a review) . Unfortunately, 
the errors are so large that it would be premature 
to infer from this a scale dependence of G. Even 
worse, phenomenologically there is a proportion- 
ality constant between the galaxy-galaxy corre- 
lation function and the actual mass correlation 
function, the so- called biasing factor, which is 
supposed to be scale dependent and could mimic 
a variable gravitational constant. However, fu- 
ture sky surveys might be able to constrain more 
strongly the relation ( JTl] ) by measuring peculiar 
velocities with better accuracy at larger distances 
and it might then be possible to extract the scale- 
dependence of G. 

5. CONCLUSIONS 

We have seen that the running of the gravi- 
tational coupling is compatible with the obser- 
vational fact that the rotation curves of galaxies 
are constant provided some amount of baryonic 
dark matter is allowed, actually about 44% less 
than what is required for a constant G. This 
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could also explain why we see less microlensing 
events towards the halo than in the direction of 
the bulge of our galaxy. Failure in reproducing 
the predicted distribution of baryonic dark mat- 
ter would signal either that the approach adopted 
here is unsuitable or that the fit (||) is inadequate. 
We have looked for possible bounds on variations 
of G with scale from primordial nucleosynthesis, 
variations in the period of binary pulsars, macro- 
scopic gravitational lcnsing and deviations in the 
peculiar velocity flows. Unfortunately, as obser- 
vational errors tend to increase with the scale 
probed, we cannot yet seriuosly constrain an in- 
crease of G with scale, as proposed by the asymp- 
totically free theories of gravity. 
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